Interface ordering and phase competition in a model Mott-insulator— band-insulator 

heterostructure 
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The phase diagram of model Mott-insulator-band-insulator heterostructures is studied using the 
semiclassical approximation to the dynamical-mean-field method as a function of thickness, coupling 
constant, and charge confinement. An interface-stabilized ferromagnetic phase is found, allow the 
study of its competition and possible coexistence with the antiferromagnetic order characteristic of 
the bulk Mott insulator. 
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Fabrication and investigation of heterostructures in- 
volving correlated-electron materials^ are an important 
direction in material science. Understanding of the elec- 
tronic properties near the interfaces and surfaces is not 
only of scientific interest but is also indispensable to 
realize electronic devices exploiting the unique proper- 
ties of correlated-electron materials. A variety of het- 
erostructures have been fabricated and studied including 
high-T c cuprates^^ Mott-insulator and band-insulator 
heterostructure^ and superlattices of transition-metal 
oxides . 6 i 7 i 8 Interestingly, the heterostructures comprising 
of a Mott-insulator and a band-insulator were reported 
to show metallic behavior, 5 

A fundamental question raised by those studies is 
"what electronic phases are realized at interfaces." For 
the vacuum-bulk interface (i.e., the surface), Potthoff 
and Nolting)* Schwieger et ai.p* and Liebschii have ar- 
gued that the reduced coordination may enhance cor- 
relation effects. The enhanced correlations could pre- 
sumably induce surface magnetic orde ring, alt hough this 
possibility was not discussed in Refs. liilfll ll Matzdorf 
et al. proposed that ferromagnetic ordering is stabilized 
at the surface of two-dimensional ruthenates by a lattice 
distortion, 12 " but this is not yet observed. Surface ferro- 
magnetism had been also discussed in a mean field treat- 
ment of the Hubbard model by Potthoff and Nolting. 13 
Similarly, the effect of bulk strain on the magnetic order- 
ing in perovskite manganites was discussed by Fang et 
alii 

All these studies dealt with systems in which charge 
densities remain unchanged from the bulk values, and 
physics arising from the modulation of charge density 
was not addressed. A crucial aspect of the recently fab- 
ricated heterostructures is charge inhomogeneity, caused 
by the spreading of electrons from one region to another. 
Strongly correlated materials typically possess interest- 
ing density-dependent phase diagrams; raising the possi- 
bility of interesting phase behavior at interfaces. In this 
paper we use the semiclassical approximation (SCA)A£ to 
the dynamical-mean-field methodic to explore the phase 
behavior of a model Mott-insulator-band-insulator het- 
erostructure. The SCA is computationally inexpensive 
and a good representation of phase diagrams and tran- 
sition temperatures in several models, allowing us to in- 



vestigate the phase behavior for a wide range of parame- 
ters, and in particular to access the T > regime which 
Hartree-Fock and related approximations fail to repre- 
sent adequately^ We observe antiferromagnetic order- 
ing in regions with charge density ~ 1 characteristic of 
a bulk Mott insulator, while ferromagnetic ordering is 
found to be a surface effect supported by an intermediate 
charge density and a strong coupling. However, we have 
found that, despite the successes noted in previous work, 
the SCA overestimates ferromagnetism on the lattice we 
study here. Therefore, our results should be regarded 
as qualitative explanations of the type of phase behavior 
which may occur rather than as quantitative statements 
about the Hubbard-modcl phase diagram. 

We study the model heterostructure introduced in 
Ref. 0. The Hamiltonian is a simplified representation 
of conduction bands of the systems studied in Ref. with 
the orbital degeneracy neglected. We consider [001] het- 
erostructures formed by varying the A-site of a ABO3 
perovskite lattice. The electrons of interest reside on the 
-B-site ions, which form a simple cubic lattice with sites 
labeled by i as fi — (xi,yi,Zi) — a(ni,rrii,li) with the 
lattice constant a set to unity. We assume each i3-site 
has a single orbital; electrons hop between nearest neigh- 
bor sites with the transfer t. The electrons interact via a 
on-site interaction U and a long-ranged Coulomb repul- 
sion. The heterostructure is defined by n planes of charge 
+1 counterions placed on the A' sublattice of A-site ions 
at positions rf = a(nj + 1/2, mj + 1/2, lj + 1/2), with 
— 00 < rij, rrij < 00 and lj = 1, . . . , n. (Coordinate z will 
be shifted such that the center of the heterostructure, A' 
sublattice, comes to z = 0.) Charge neutrality requires 
that the areal density of electrons is n. The resulting 
Hamiltonian is H = H band + H rnt + H Co ui with 
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Note that [/ 7^ on all sites. A dimensionless measure 
of the strength of the long-ranged Coulomb interaction 
is E c — e 2 /(eat) with the dielectric constant s. In most 
of our analysis, we choose E c = 0.8. This corresponds 
to t ~ 0.3 eV, b~4A, and e — 15, which describe the 
system studied in Ref. Oj. The charge profile is found not 
to depend in an important way on e, but the stability 
of magnetic orderings does because this is sensitive to 
the details of the charge density distribution as discussed 
later. 

The basic object of our study is the electron Green's 
function. In general, this is given by 

Gcr(r, r'; u) = [w + fi- H band - H Cm d - S CT (r, r'; 

(4) 

with the chemical potential fi and the electron self-energy 
S. We consider [001] heterostructures with either in- 
plane translational invariance or Ag-sub-lattice antiferro- 
magnetism. The Green's function and self-energy are 
therefore functions of the variables (z, rj, z', r]', fcii) where 
r\ and rj'(— 1, . . . , N s ) label the sublattice in layers z and 
z', respectively, and k\\ is a momentum in the (reduced) 
Brillouin zone. As in Ref. 0, we approximate the self- 
energy as the sum of a static Hartree term arising 
from the long-ranged part of the Coulomb interaction and 
a dynamical part (ui) arising from local fluctuations. 
Generalizing the inhomogeneous case of the dynamical- 
mean- field theory (DMFT)jifii we assume that the self- 
energy is only dependent on layer z and sublattice 77. 
Thus, the dynamical part of the self-energy is written as 

=*£?(z,»j;w). (5) 

The ^-dependent self-energy is determined from the so- 
lution of a quantum impurity modeU^ with the mean- 
field function fixed by the self-consistency condition 

G^(z, m u)^N s G a (z, V ,z, V ,k f cu). (6) 

We shall use the DMFT, which is a mean-field ap- 
proximation, to calculate magnetic phase diagrams of ef- 
fectively two dimensional systems, so some mention of 
fluctuation effects is needed. In general, the calculated 
transition temperatures are to be understood as crossover 
scales below which the magnetic correlation length £ 
grows rapidly £ ~ exp(27rp s /T) with spin stiffness 
discussed in more detail below. True long-ranged order 
may be induced at T > by an Ising anisotropy or by 
coupling in a system made of repeated heterostructure 
units. In any event, the rapid growth of the correla- 
tions for T < p s means that the properties are effectively 
those of an ordered state. We will find that the low-T p s 
is sufficiently large relative to the calculated transition 
temperature that fluctuation effects are not crucial. 

In general for the heterostructure with L layers with 
N s sublattices, one must solve L x Ns independent im- 
purity models. Due to the self-consistency condition [cf. 



Eq. @] and to compute the charge density n a (z,r}) = 
- J ^/(w)ImG™ p (z, 77; ui) with / the Fermi distribution 
function, it is required to invert the (L x Ns) 2 Green's 
function matrix at each momenta and frequency. This 
time consuming numerics restricts the size of the unit 
cell. In this study, we consider the commensurate mag- 
netic states with up to two sublattices, N s = 1 and 2, 
on each layer and with the charge density independent of 
the sublattices, i.e., paramagnetic (PM), ferromagnetic 
(FM) states, and (layer-) antiferromagnetic (AF) state 
where antiferromagnetic (FM) planes with moment al- 
ternating from plane to plane. Note that the AF state 
extrapolates to the bulk AF state with the magnetic vec- 
tor q = (7r,7r,7r) at n — > 00. By symmetry, the number 
of quantum impurity models one must solve is reduced 
to L. However, solution of the impurity models is a time 
consuming task, and an inexpensive solver is required. In 
Ref. 0, to study the evolution of the low-energy quasi- 
particle band and high-energy Hubbard bands as a func- 
tion of position, we applied two-site DMFT— which is 
a simplified version of exact-diagonalization method. At 
T = 0, this method is known to give reasonable result for 
Mott metal-insulator and magnetic transitions. However, 
small number of bath orbitals is known to be insufficient 
to describe the thermodynamics correctly^ 

For the investigations presented here, we use the semi- 
classical approximation, which is computationally inex- 
pensive and has been found to be reasonably accurate 
for phase boundaries and excitation spectra of several 
models, including the half-filled Hubbard model and, at 
all fillings, for the d — 3 and d — 00 face-centered-cubic 
(FCC) lattices^ We note however that, perhaps because 
it does not properly include quasiparticle coherence, the 
method overemphasizes ferromagnetism at intermediate 
density, giving for the FCC lattices transition tempera- 
tures ~ 50 % higher than those found by quantum Monte- 
Carloji^ and for cubic lattices, finding ferromagnetism at 
ntot — 0.5 and moderately large U (of order the criti- 
cal value for the Mott transition) when other methods 21 
suggest ferromagnetism is confined to very large U and 
ntot near 1. Thus, the application of the SCA to Hub- 
bard heterostructure allows a convenient exploration of 
the interplay between different phases, but probably does 
not provide a quantitatively reliable picture of the phase 
diagram. 

First, we investigate the magnetic behavior at finite 
temperature. The upper panel of Fig.^shows our calcu- 
lated phase diagram in the interaction-temperature plane 
for heterostructures with various thicknesses for charge 
binding parameter E c = 0.8. The one-layer heterostruc- 
ture is PM at weak to moderate interactions, and FM 
at strong interactions. The two- and three-layer het- 
erostructures are AF at weak to intermediate interaction, 
and become FM at stronger interactions with almost the 
same Tc for n — 2 and 3. Antiferromagnetic Neel tem- 
perature T/v is found to be strongly dependent on the 
layer thickness; it increases with the increase of layer 
thickness. These T/v's are substantially reduced from 
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FIG. 1: Upper panel: Magnetic transition temperatures of 
heterostructures with various thicknesses n indicated as func- 
tions of interaction strength. E c = 0.8. Filled symbols: Fer- 
romagnetic Curie temperature Tc, open symbols: Antiferro- 
magnetic Neel temperature Tn- Note that, where both phases 
are locally stable, the phase with higher Tc(Xjv) has the lower 
free energy. Lower panel: Magnetic transition temperatures 
and charge density n to t at z = (light solid line) and ±1 
(light broken line) at T/t — 0.1 as functions of the parameter 
Ec for 2-layer heterostructure with U/t — 20. Counterions 
are placed at z — ±0.5. 



the the bulk values; maximum value Tff ax /t ~ 0.47 at 
U/t ~ 10, whereas Curie temperature Tc's are almost the 
bulk 2d values at n to t ~ 0.5. The strong dependence of 
Tjv on thickness may be understood from the bulk phase 
diagram; antiferromagnetism is stabilized only very near 
to half-filling, and in the thinner heterostructures the 
charge-spreading effect reduces the density too much. 
This physics is seem from a different point of view in 
the lower panel of Fig. ^ Hartree-Fock studies of this 
and related models find a layer-AF phase. This phase is 
not found in our DMFT analysis. 

The lower panel of Fig. ^ presents a detailed study 
of the n = 2 heterostructure showing how changes in the 
charge confinement parameter E c affect the physics. The 
filled and open points (left-hand axis) show the variation 
of the Curie and Neel temperatures, respectively. The 
light solid and light broken lines (right-hand axis) show 
the variation of charge density on the central and next to 
central layers, respectively. It is seen that the AF order- 
ing is rapidly destabilized with the decrease of E c (and 
concomitant decrease of charge density), and T/v is seen 
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FIG. 2: Magnetization density of 4-layer heterostructure. 
Counterions are placed at z — ±0.5, ±1.5. E c — 0.8. Up- 
per panel: Magnetization m in a FM state for U/t = 24 and 
T/t = 0.1. Lower panel: In-plane staggered magnetization 
m s tagg in an AF state for U/t = 6 and T/t = 0.05. Filled 
(open) circles are the results by DMFT (HF). For compari- 
son, charge density ntot computed by DMFT is also shown 
in the lower panel (filled squares). Note that the staggered 
magnetization in the outer layers (\z\ > 2), and the outermost 
layers (\z\ = 1) have the same sign. 



to be correlated to that of the charge density at z = 0. 
On the contrary, Tq has a weak variation with E c . This 
indicates that the FM ordering is favored by the interme- 
diate charge density as discussed in the bulk single-band 
Hubbard model; 22 at large E c , the magnetization is large 
on the outer layer and small in the inner layer, at small 
E c , the situation is reversed. 

These arguments are confirmed by our calculations of 
the spatial variation of the magnetization density, re- 
ported in Fig. [21 which shows numerical results for a 4- 
layer heterostructure with counterions at z = ±0.5 and 
±1.5, and E c = 0.8. The upper panel of Fig. [2] shows the 
magnetization in the FM state. In DMFT (filled circles), 
only the layers near the interfaces (\z\ ~ 1-2) have large 
polarization and inner layers in the heterostructure have 
small moments. This explains the weak n-dependencc 
of Tc of thick heterostructures (see the upper panel of 
Fig.[IJ. In HF (open circles), all layers in the heterostruc- 
ture are highly polarized. In contrast in an AF state, 
result of the in-plane staggered magnetization by DMFT 
and HF agree well as shown in the lower panel of Fig. [21 
For comparison, the total charge density is also plotted 
(filled squares). The in-plane staggered magnetization is 
large only at inner layers where the charge density is close 
to 1. Note that the staggered magnetization in the outer 
layers (\z\ > 2) has the same sign as in the outermost 
layers (\z\ = 1) indicating that the outer layers are not 
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FIG. 3: Layer- and sublattice-resolved spectral functions as 
functions of real frequency ui for 4-layer heterostructure. E c — 
0.8. Upper panel: Ferromagnetic state at U/t = 24 and T/t — 
0.1. Lower panel: Antiferromagnetic state at U/t = 6 and 
T/t = 0.05. Sublattice 77 = 1. Solid (broken) lines are for up 
(down) spin electrons. For sublattice rj = 2 in AF state, up 
and down electrons are interchanged. 



intrinsically magnetic. 

We now return to the issue of fluctuation effects. We 
study a model which is two dimensional and order param- 
eters with a continuous spin rotation symmetry; thus, at 
T > 0, rather the calculated Tc means a crossover to a 
low-T "almost ordered" phase characterized by a expo- 
nentially growing correlation length £ ~ exp(27rp s /T) 19 
where the key stiffness p s is given by the second deriva- 
tive of the free energy with respect to the order parameter 
orientation. We have computed p s for a single ferromag- 
netic plane with charge density 0.5 finding p s ~ i/8 We 
observe that in all cases several planes exhibit the rele- 
vant order, so that the total stiffness is larger by a fac- 
tor 2-3. Thus we see that the stiffness is large enough 
relative to the DMFT Jp^'s that the purely two dimen- 
sional fluctuation effects are of minor importance, leading 
to £ > exp(4-5) for T < T c , N /2. 

Spin distributions presented in Fig. [21 can be under- 
stood from the single-particle spectral functions. In 
Fig. |3 are presented the DMFT results for the layer- 
and sublattice-resolved spectral functions A a (z,n;uj) = 
-±ImG i Z np (z > r);uj + i0+) for the FM (upper panel) and 
the AF (lower panel) states of 4-layer heterostructure 
with the same parameters as in Fig. These quantities 
can in principle be measured by spin-dependent photoc- 
mission or scanning tunneling microscopy. As noticed 
in Ref. ILsl spectral function outside of the heterostruc- 
ture (\z\ ^> 2) is essentially identical to that of the free 



tight-binding model Hb an d, and electron density is neg- 
ligibly small. With approaching the interfaces (\z\ = 2), 
the spectral function shifts downwards and begins to 
broaden. In the FM case, magnetic ordering is possible 
only near the interface (\z\ ~ 1-2) carrying the intermedi- 
ate charge density. Inside the heterostructure (\z\ <C 2), 
clear Hubbard gap exists due to the large U and uniform 
polarization is hard to achieve. On the contrary, high 
charge density is necessary to keep the staggered magne- 
tization in the AF case as seen as a difference between 
up and down spectra in the lower panel of Fig. [21 

So far, we have discussed the competition between fer- 
romagnetic and antiferromagnetic phases and spatial dis- 
tributions of charge and magnetic densities establishing 
that ferromagnetism is a surface effect. Finally, we dis- 
cuss the interplay between phases. We note that the 
charge density or magnetization density typically varies 
over ~ 3 unit cell range, so for thin heterostructures and 
moderate charge confinement energies, the entire het- 
erostructure exhibits a single phase, controlled by the 
instability exhibiting the highest transition temperature 
(cf. Fig. However, for thicker heterostructure or 
stronger confinement we believe that one can observe an 
ordered state involving an antiferromagnetic center with 
a ferrimagnetic "skin." A hint of this behavior can be 
observed in Fig. [21 Consider the central layer; this has 
an occupancy n to t ~ 0.9, so from Fig. ^ and rescaling 
account for the different U, we would expect Tjy ~ O.lt. 
Now this layer is not isolated; hoppings to the layers of 
z = ±1 leads to a effective polarizing field of the or- 
der of 2(1 — n to t)t m I ntot where the first factor is the 
number of layers, the second is the hopping amplitude 
renormalized by strong correlations, and the third fac- 
tor is the relative spin polarization; putting these factors 
together gives a polarizing field of about 0.15t, approx- 
imately equal to the AF coupling. Thus for the central 
layer of this heterostructure ferromagnetic and antifer- 
romagnetic tendencies are very closely balanced, but for 
thicker systems (not at present computationally accessi- 
ble) or perhaps for stronger charge confinement, an anti- 
ferromagnetic center will occur. 

To summarize, we have presented a semiclassical 
DMFT study of magnetic phase behavior of a model 
Mott-insulator-band-insulator heterostructure in which 
the behavior is controlled by the spreading of the elec- 
tronic charge out of the confinement region. Magnetic 
phase diagram is investigated as a function of layer thick- 
ness, temperature, and interaction strength. Ferromag- 
netic ordering is found to be a surface effect stabilized at 
an interface region with moderate charge density, while 
antiferromanetic ordering is found at a region with high 
density ~ 1 characteristic of the bulk Mott insulator, 
and Neel temperature is sensitive to the layer thickness 
and charge confinement energy. These magnetic order- 
ings may coexist in very thick heterostructure exhibiting 
a ferromagnetic "skin" and an antiferromagnetic "core." 
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